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Thermally activated magnetization reversal in elongated particles is studied within a model that allows for spatially nonuniform 
magnetization configurations along the particle. An external field antiparallel to the existing magnetization is shown to give rise 
to an energy barrier which represents a spatially localized deviation from the initial uniform magnetization configuration. For 
sufficiently elongated particles thermal fiuctuations thus substantially lower the coercivity compared to the previous theories 
by Neel and Brown which assume a spatially uniform magnetization. The magnetization reversal rate is calculated using a 
functional Fokker-Planck description of the stochastic magnetization dynamics. Analytical results are obtained in the limits of 
small fields and fields close to the anisotropy field. In the former case the hard-axis anisotropy becomes effectively strong and 
the magnetization reversal rate is shown to reduce to the nucleation rate of kink-antikink pairs in the overdamped sine-Gordon 
model. The present theory therefore includes the nucleation theory of the double sine-Gordon chain as a special case. 

PACS numbers: 05.40.+j, 75.10.Hk, 75.60.Ch, 75.60.Ej, 75.60.Jp, 85.70.Li, 87.40.-fw, 91.25.Ng 



I. INTRODUCTION 

The magnetization in a uniformly magnetized sample 
is usually stabilized by an easy-axis anisotropy of crys- 
talline or demagnetizing origin. To reach a state of zero 
net magnetization one has to apply an external field in 
the reversed direction, the so called coercive field. In 
macroscopic samples of high purity such as Yttrium Iron 
Garnet (YIG)l!l, this field can be less than 10"^ Qc. This 
low coercivity is commonly attributed to the existence of 
residual domains of reverse magnetization in the original 
uniformly magnetized state. The measured coercivity is 
then associated with the depinning and motion of the 
corresponding domain walls. 

This situation is strikingly different for microscopic 
single-domain particles where no such residual domains 
exist. The coercivity can reach values of more than 1000 
Oe since the state of reversed magnetization has first to 
be nucleated. Consequently such particles exhibit an ex- 
tremely high long-term stability of the magnetization. 
This fact renders them suitable for information storage 
in recording media and as constituents of rocks they pre- 
serve the value of the local magnetic field as the temper- 
ature has dropped below the blocking temperature of the 
particle. With decreasing sample size, however, the effect 
of thermal fluctuations becomes increasingly important. 
For particle sizes of a few nanometers and at room tem- 
perature the magnetization fluctuates randomly over the 
anisotropy barrier and a superparamagnetic state results 
with vanishing average magnetic moment. 

This paper concentrates on particles whose size is 
above the superparamagnetic limit but which are still 
small enough that the coercivity is affected by thermal 
fiuctuations. iThc only ab-initio theory has been devel- 
oped by NeelEl and BrownH and it is based on the as- 
sumption that the magnetization distribution is uniform 
throughout the sample. Consequently the energy bar- 



rier is proportional to the volume and the Arrhenius fac- 
tor leads to an exponential suppression of thermal ef- 
fects with the particle volume. This picture is indeed 
adequate for small particles of approximately spherical 
shape. However, for sufficiently elongated particles a 
magnetization reversal via a rigid rotation of the mag- 
netization becomes energetically unfavorable. It will be 
more advantageous to form a spatially localized excur- 
sion from the metastable state since the additional cost 
of exchange energy due to the spatial nonuniformity is 
by far outweighed by the gain of anisotropy energy by 
keeping the deviation localized. 

It is the purpose of this paper to formulate an ab ini- 
tio theory of this effect and to show that for sufficiently 
elongated particles a spatially nonuniform barrier yields 
a much lower coercivity than previous theories. A short 
account of the results, of the present paper has already 
been given elsewhereQ. We shall start from a classical 
one-dimensional model energy density which takes into 
account the exchange interaction between the magnetic 
moments along the particle. In addition, the energy den- 
sity contains hard- and easy axis anisotropics as well 
as the coupling to an external field. The anisotropics 
may contain contributions of both, shape and crystalline 
anisotropics. The barrier energy is then shown to be in- 
dependent of the hard-axis anisotropy, and it is propor- 
tional to the domain wall energy and the sample cross 
sectional area. Consequently, the barrier energy is inde- 
pendent of the particle length for sufficiently elongated 
particles. 

In order to induce magnetization reversal, thermal fluc- 
tuations have to form a "nucleus" of critical size with 
the property that smaller deformations fall back to the 
metastable state whereas larger deformations grow with 
energy gain until the magnetization is reversed. There- 
fore the nucleus represents an unstable structure with 
exactly one unstable mode. An analytical expression for 
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this structure has been obtainedoH and the spin wave 
excitations of |tlie nucleus have been investigated in a 
previous papercl (henceforth referred to as I). For exter- 
nal fields close to the anisotropy field which renders an 
individual magnetic moment unstable, the nucleus rep- 
resents only a small deviation from the metastable state. 
For small fields, the nucleus consists of two well separated 
domain walls enclosing an already reversed domain. 

The present approachjelies on methods that have been 
used in the descriptionQ of the dynamics of first order 
phase transitions. This method has been applied for 
the description of the decay of a supercurrentS in a thin 
wire or the propagationO'Ed of dislocations. In contrast to 
these applications we consider here the regime of moder- 
ate damping since damping in magnetic systems is very 
small. The rate is shown to be the product of a pref- 
actor depending on the external field and temperature 
T, and the Arrhenius factor exp{— ylfs/fc^T} which in- 
volves the barrier energy A£s with A the sample cross 
sectional area (fcs is the Boltzmann constant). 

For the evaluation of the prefactor we shall employ two 
different approaches. One is the Jacobi method which 
relies on the explicit knowledge of a zero-energy (Gold- 
stone) mode. Therefore it can only be applied for easy- 
plane fluctuations but it cannot be used for out of easy- 
plane excitations due to the existence of a mass gap. The 
second method makes use of the scattering phase shifts 
of spin waves around the nucleus. This latter approach 
reveals some considerable subtleties which do not seem to 
have been noted previously. First, for even-parity wave- 
functions the ID versionlil of Levinson's theorem alters 
the usual expression for the density of states. Second, 
the number of bound states of the fluctuation operators 
is not conserved under small and large nucleus approx- 
imations. This fact raises doubts about the commonly 
employed approach of performing functional integrals of 
the free energy after already having performed the limits 
of small or large nuclei. By a careful investigation it is 
shown that these two subtleties conspire in such a way 
that this interchange of limits is indeed legitimate. 

There are only sparse treatments of magnetization 
reversal in the literature. However, this field is 
closely related to macroscopic quantum tunneling of the 
(sublattice-) magnetization in small (anti-)fcrromagnetic 
grains, a subject that has attracted much interest re- 
cently. It thus appears to be useful to relate some impor- 
tant papers that contributed to the development of these 
fields. 

Early work pa nucleation theory culminated in the cel- 
ebrated papenlj by Kramers who calculated the escape 
rate due to thermal activation out of a metastable state 
in the limits of low as well as moderate to large fric- 
tion. He showed that the rate is given by a prefactor 
times the Arrhenius term. Despite the fact that his work 
was restricted to one degree of freedom, his method of 
the evaluation of the prefactor turned out to be so pow- 
erful that its spirit still underlies much more complex 
applications. An extension to an arbitrary number of de- 



grees of freedoms in the large friction limit, was due to 
Brinkman, Landauer-Swanson and LangenU. The, case 
of moderate friction has been considered by Langeiul who 
also pointed out that the nucleation rate may be inter- 
preted as the analytic continuation of the partition func- 
tion. This idea is closely rjelated to the subsequently de- 
veloped instanton conceplfij in Euclidean quantum field 
theories. Kramers's theory and its extensions have pe.- 
cently been reviewed by Hanggi, Talkner and Borkovedij. 

The first application of Kramers's theory to magnetic 
systems has been made by BrownB who investigated ther- 
mally activated uniform magnetization reversal in small 
ferromagnetic particles to explain superparamagnetism. 
He set up the Fokker-Planck equation for the stochastic 
dynamics of the masrietization and thus related Neel's 
earlier considerationa3 on reversal rates with the general 
framework of statistical mechanics. For axial symmetry 
of the anisotropy he obtained nucleation rates from the 
lowest nonzero eigenvalue of the Fokker-Planck equation 
in the limit of low barriers. For a high barrier he used 
Kramers's procedure to evaluate the rate constant. Later 
the lowest positive eigenvalue was investigated numeri- 
cally for all intermediate values between low and large 
barriers by AharoniliH. Eisenstein and Aharonilia inves- 
tigated the competition of the uniform mode and the 
nonuniform curling mode as possible candidates of crit- 
ical nuclei for different particle radii. However, the nu- 
cleation rates via thcunonuniform mode were calculated 
using Brown's theorytl for spatially uniform nucleation. 

Subsequently the issue of magnetization reversal rates 
was not addressed for many years. A renewal of inter- 
est then arose from a quantum mechEpical point of view. 
Based on path integral and WKBt3 techniques, first 
investigations showed that a single spin in an anisotropic 
field behaves similar to a particle with inertia and tun- 
nels between difierent anisotropy minima. It has then 
been suggestedEil that in small ferromagnetic particles 
macroscopic quantum tunneling might occur. The effect 
of dissipation due to magnetoelastic coupling has been 
discussed by-parg and KimEHl. In the context of recent 
experimentaH3 these apprQa|Ches have been reexamined 
and it has been predictedEj that quantum tunneling is 
suppressed for half integer spins as a consequence of a 
previously neglected Wess-Zumino term in the quantum 
spin action and the destructive interference of instanton 
and antiinstanton paths. 

While all these approaches dealt with tunneling via 
spatially uniform structures, tunneling via spatially 
nonuniform (bubble) structures in two dimensions was 
investigated in the liraiti of external magnetic fields close 
to the anisotropy fieldciand for very small fields in the 
thin wall approximations. In the latter case the nucle- 
ating structure is a large cylindrical domain of reversed 
magnetization delimited by a Bloch wall. Various as- 
pects of quantum tunneling with emphasis on tunneling 
of Bloch walls have recently been reviewed by Stamp, 
Chudnovsky and BarbaracJ. 

Surprisingly, the conceptually much simpler classical 
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problem of thermal nucleation remained untouched un- 
til recently. Klik and Guntheicj calculated the nucle- 
ation rate for nucleation via uniform structures for cubic 
symmetry. In contrast to earlier investigations they also 
calculated nucleation rates for a weakly damped SKStem. 
Nuclei of curling symmetry in an infinite csdindei£3 and 
a nucleation center of spherical symmetrytil have been 
investigated recently. 

From this review, there emerges clearly the need of 
an ab initio theory for magnetization reversal rates via 
spatially nonuniform structures. The present work is or- 
ganized as follows: 

In section II some results of paper I are reviewed which 
are relevant to the present work. In section III a func- 
tional Fokker-Planck equation is constructed which de- 
scribes the stochastic magnetization dynamics near the 
nucleus and the corresponding nucleation rate is derived. 
It is shown that the result has the same general struc- 
ture as that of Ref . ^b) . The prefactor separates into a 
term describing the dynamical decay of the nucleus and 
in a term arising from the Gaussian fluctuations around 
the nucleus. The unstable mode enters in such a way 
as if it represented a stable mode. The details of the 
calculations are presented in an appendix. The explicit 
evaluation of the prefactors in various limits is then the 
subject of the remaining sections. In section IV we eval- 
uate the statistical part of the prefactor analytically in 
the limit of small and large nuclei as well as for large and 
small values of the hard-axis anisotropy. For small nuclei 
and if the hard-axis anisotropy is much larger than the 
easy-axis anisotropy, the out of easy-plane fluctuations 
do not contribute at all. In section V the nucleation rate 
is evaluated in the overdamped limit. The rate in the 
moderately damped regime and the decay frequency of 
the nucleus are investigated in section VI. In section VII 
the results of the previous sections arc used to calculate 
the creation rate of kink-antikink pairs in the double sine- 
Gordon system. It is shown that this rate reproduces the 
magnetization reversal rate in the limit of large hard- 
axis anisotropy or external fields close to the (easy-axis) 
anisotropy field. In section VIII experimental implica- 
tions are discussed. For a particle of lOOA diameter and 
an aspect ratio 15 : 1, the present theory is shown to yield 
a coercivity reduction from the anisotropy field that is 
twice as large as that of the Neel-Brown theory. Finally 
the applicability range of the present theory is discussed 
since it is known that in the underdamped limit the rate 
is governed by a diffusion in energy rather than in con- 
figuration space. 



II. MODEL, NUCLEUS AND FLUCTUATIONS 

In this section we present the model and review some 
important results of paper I. The ferromagnet is de- 
scribed within a classical continuum model, the magneti- 
zation being represented by a vector M of constant mag- 



nitude Mq. We focus on an effectively one-dimensional 
situation where the magnetization only depends on one 
coordinate, i.e. M — M(a;,i). The energy per unit area 
is given by 



L/2 c . 

-L/2 I ^^-^O 



(2.1) 



where dx = d/dx and L is the finite sample length in x- 
direction. Ultimately, we shall be interested in the limit 
L ^ oo. The first term in (2.1) is the classical coun- 



terpart of exchange energy and A is an exchange con- 
stant. The second term defines an easy-axis along the 
x-direction. The third term is a hard-axis anisotropy 
which favors the magnetization to lie in the a;?;-plane. 
Ke > and Kh > are easy- and hard-axis anisotropy 
constants respectively. The degeneracy between the two 
anisotropy minima along the a;-axis is lifted by an ex- 
ternal magnetic field -ffext along the positive x-axis. The 
energy (2.1) describes magnetization configurations in an 
elongated particle of diameter smallero or comparable to 
the minimal length scale in the system A/ K^^y^, where 
^^max is the larger of the anisotropy constants K^, Kh- 

Note that the energy (2.1) can be used to describe 
three distinct anisotropy configurations in elongated par- 
ticles. The first, most common case is an easy- axis along 
the particle axis which may be caused by both demag- 
netizing (shape) and crystalline anisotropy (cf. Fig |l| 
a). E.g. for an infinite cylinder with an easy-axis along 
the sample one has Ke — t^Mq + -fiTe, cryst where the first 
term is due to the shape anisotropy and the second term 
due to crystalline anisotropy. The hard-axis anisotropy 
may arise either from an additional crystalline easy-axis 
that is misaligned with the particle axis or from an el- 
liptic sample cross section. The second case is an elon- 
gated particle of a material of high crystalline anisotropy 
(i^e^cryst > 2ttMq) with both easy- and hard-axis per- 
pendicular to the long- axis of the sample (cf. Fig |l|b). 
The third case refers to a thin slab with an easy-axis 
anisotropy in the film plane (cf. Fig|^ c). 

In the following we focus on a situation as in Fig. la). 
The results for configurations shown in Figs. 1 b),c) 
are simply obtained by substituting y and z for the x- 
dependence of the magnetization. The components in 
internal (spin) space remain unchanged and the spheri- 
cal coordinates are always defined in the same way with 
respect to the coordinate axes. 

The dissipative dynamics of the magnetization is as- 
sumed to obey the Landau-Lifshitz-Gilbert equations 
(see e.g. Ref. |3^ ): 

dtM = -7M X Heff + — M X dtM. (2.2) 

Mq 

where 7 > denotes the gyromagnetic ratio, a > 
is the dimensionless damping constant and dt = d/dt. 
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The first term on the r.h.s. of (2.2) describes the preces- 
sion of the magnetization in the effective magnetic field 
Hefi = — (Jf/^M ((5/( ?M denotes a functional derivative). 
The second term in (2.2) is a viscous damping term and 
accounts for the relaxation of the magnetization into the 
direction of the effective magnetic field. This term is phc- 
nomenological in nature. It describes damping processes 
which conserve the magnitude of the magnetization at 
every space point. It is convenient to rewrite (2.2) as 
follows 



(l + a^)5tM = -7MxHeff 



q;7 

Mq 



M X [M X Hoff] 



(2.3) 



This equation i s ob tained by evaluating the cross prod- 
uct of M with ( ^.2|) . Eq. (2^) is similar to the damping 
term origi nally proposed by Landau and Lifshitz. How- 
ever, Eq. (2.3) contains the damping parameter a such 
that the motion is slowed down for large a while the orig- 
inal equation of Landau-Lifshitz exhibits an unphysical 
accelerat ion of the motion for large damping parameters. 

Since (2.3) conserves the magnitude of the magnetiza- 
tion it is appropriate to introduce spherical coordinates 
according to M/A/q = (sin cos 0, sin sin 0, cos 6*) . We 
use dimensionless units defined by 




[t] =il + a') 



[£] = 



(2.4) 



\J AjKf. is the Bloch wall width, is the preces- 

sion frequency in the anisotropy field. To simplify nota- 
tion, an additional factor l-\-a^ has been absorbed in the 
time scale. 2\/~AKl is half the energy of a static 7r-Bloch 
wall. In dimensionless units and spherical coordinates 



the energy per area (2.1) takes the following form 



[sin^ cos^ 



J-L/2 {2 

^ - 1] + - — cos^ e - h[sm 0co8(l)+ 1] 



In ( |2.5| ), the dimensionless anisotropy ratio 



and the reduced external field 

n = . 

2K, ' 



(2.6) 



(2.7) 



have been introduced. Using ( |2.4|) a nd spherical coordi- 
nates, the equations of motion ( ^.3^ t ake the following 
form (see also the appendix of Ref. |33|) 



sin 9 dt4 



19 



a 



1 se: 

sin 9 56 



1 6£ 5£_ 

sk^du" "19 



(2.8) 



The first terms on the r.h.s. describe the precession in the 
effective magnetic field, whereas the terms proportional 
to a are damping terms. 

Spatially uniform static solutions of ( |2.3| ) lie in the 
easy plane and are given by {(po,9o) = (0, 7r/2) and 
{<l}m,dm) — (tt, 7r/2), the latter being stable only for 
h < 1. The state (0o,6'o) is completely aligned with 
the external field and thus represents the state of lowest 
energy. The configuration {9m, (pm) is oriented antipar- 
allel to the external field and its energy per volume ex- 
ceeds that of the ground state (0Oi^o) by 2h. Therefore 
) is a metastable state for ft, < 1. 

At finite temperatures, the magnetization exhibits fluc- 
tuations around the metastable state until it eventually 
overcomes a barrier for magnetization reversal. For large 
sample lengths L, a magnetization reversal via a uniform 
rotation of the magnetization is highly unlikely since it 
would require an energy proportional to L. Instead, the 
system will establish magnetization reversal by forming 
a spatially localized deviation from the metastable state. 
There is a well defined "nucleus" of critical size with the 
property that deformations of smaller size fall back to 
the metastable state, whereas larger deformations grow 
with energy gain until the whole system is in the ground 
state parallel to the external field. In paper I it has been 
shown that the magnetization configuration (cf. (1,3.9)) 
defined by 



0± cosh ^ 

tan-^ = ± 9s^tt2, 

2 sinhi? ' ^ / ' 



with 



sech R ^ h. 



S = coth R, 



(2.9) 



(2.10) 



exhibits exactly one unstable mode an d th us represents 
such a nucleus of critical size. Eq. ( |2.9|) is in princi- 
ple only valid for a sample of infinite length L but it is 
an excellent approximation for a sample of finite length 
(2.5^ L > 21:^ AjKe- For typographical simplicity the sub- 
script s of (5 and R has been dropped in contrast to paper 
/. denotes the arbitrary position of the nucleus along 
the particle. This degeneracy with respect to translations 
gives rise to zero energy (Goldstone) mo de. I n the follow- 
ing we shall put = 0. The structure ( |2.9| ) can also be 
viewed as a superposition of two 7r-Bloch walls centered 
at X = ±i?/(5 + x q with opposite relative sense of twist. 
For R small, (^|^) represents only a small deviation from 
the metastable state, whereas for large R it represents a 
large domain of size 2R5 delimited by an untwisted pair 
of domain walls (cf. Fig. ||). In the following we shall 
restrict ourselves to 0^ only and we shall drop the su- 
perscript. The existence of two equivalent saddle points 
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will result in a factor of two in the final expression of the 
nucleation rate. 

Out of easy-plane fluctuations p and azimuthal fluctu- 
ations around are introduced as follows 



{x,t) = 7r/2 



(2.21) 



into (2.5) we have 



{X,t) = (j)s{x) +ip{x,t), 
{x,t)=TT/2-p{x,t). 



(2.11) 

nd 



L/2 



L/2 



dxipTL'- 



'P + t; 



L/2 



L/2 



dx p H'^P p, 



Inserting ( 2.11D into the energy ( |2.5[ ) we obtain up to 2 
order in ip and p 



(2.22) 



£:(2) - P + _ 



L/2 



L/2 



dx f 



^+2 



L/2 



L/2 



where the operators 
dxpU'Pp, (2.12) 



where 



dx^ 
WP = 

dx'^ 



(2.13) 
(2.14) 



The energy per unit area of the nucleus (2.9) is given by 
f,, =4tanhi?-4i?sech^i?. (2.15) 



The characteristic width 5 is given by ( 2.1C ), and the 
potentials V± can be inferred from (1,4.13) 

V±(C, i?) = 1 - 2 sech^(^ + R)~2 sec\?{e, - R) 

±2sech(^-Ki?)sech(C-i?). (2.16) 



The eigenvalue problems of ( 2.13 ) and ( 2.14 ) are written 
as follows 

W^X7{^. R) = Et^{R)x7{x, i?), (2.17) 
H'^xfix, R) = E:p{R)x7{x, R), (2.18) 



where v denotes both, bound states and scattering s tates . 
From I, two solutions of the eigenvalue problems ( 2.17 ) 
arc known: the ground state of T-C^p 

Xo" oc (5-1 {sech(^ +R)+ sech(| - R)} , E^" - Q~\ 

(2.19) 

and the first excited state of H^"^ 



Sip 

Xi « • 



-i{sech(|-fi?)-sech(|-i?)}, El^^O. 



(2.20) 

Since x'l' is antisymmetric, there is exactly one unstable 
mode of negative energy in ip while all fluctuations in p- 
direction have positive energy since > 0. Therefore 
the untwisted domain wall pair represents a saddle point 
of the energy with exactly one unstable direction. 

The operators characterizing the modes around the 
metastable state {4>m,Qm) are obtained in an analogous 
way. Inserting 



dx'^ 



(2.23) 



reflect the spatial uniformity of the metastable state. In 
order to calculate nucleation rates by thermal activation, 
we have to examine the stochastic dynamics around the 
nucleus. 



III. STOCHASTIC MOTION AND NUCLEATION 
RATE 

The dissipative dynamics of the magnetization is gov- 
erned by the equations of motion (2.8 ). To inv estiga te the 
dyna mics near the nucleus, we insert (2.11) and (2.12) 
into (2.8) and obtain the linearized equations of motion 



dtp = -WPp - aW'^ip, 
dtp =W^Lp~ aWPp. 



(3.1) 



Notice the unusual occurrence of damping terms propor- 
tional to a in the equations of motion for both ip and p. 
The reversible part of (3.1) is of Hamiltonian structure. 
This is due to the fact that the z-component of angular 
momentum cos 9 and the azimuthal angle (p are canoni- 



cally conjugate variables. However, the signs in Eq. (3.1 



are reversed compared to a usual canonical theory since 
we are dealing with the magnetization rather than the 
angular momentum. 

The statically unstable mode has a dynamical coun- 
terpart ((p+(x, t),p+(x, t)) oc e'^+*((j9+(x), p+{x)) with 



A+ > which inserted into (3.1) obeys the coupled eigen- 
value problem: 



\+p+^-WPp+-aW'^p+, 

X+P+ = n''^ip+ - an'Pp+, 



(3.2) 



with the boundary conditions p+{±^) = p+(±-|) — 
(/j^(±-|) =p^(±-|) = 0. The linearized equations (3.1) 
of motion can also be cast in a compact form 



dti^i{x,t) = -^Mi.j'Hjtpj{x,t). 

i 

In (|3.3|) we have introduced 



(3.3) 



5 



and the dynamic matrix 



M : 



a 1 
-1 a 



which is the sum of a symplectic matrix describing the 
reversible part of the dynamics and a diagonal posi tive 
definite dissipative matrix. For the operators in (3.3) we 
have used the notation (^1,7^2) = in'"^,n'P). 

Eqns ( t3.3[) , and equivalently (^), describe the deter- 
ministic motion of the system in the vicinity of the saddle 
point. However, they are not consistent with the fluc- 
tuation dissipation theorem since they lack the stochas- 
tic forces resulting from the coupling to the heat bath. 
Without these stochastic forces, the magnetization would 
never be driven away from the initial metastable state. 
Stochastic forces can be added to the r.h.s. of (pTl|) or 



(3.4) arising from the condition that / lyipDp g^q = 1 in the 
vicinity of the metastable state. Since g^q is sharply 
peaked around the metastable state, a Gaussian approx- 
imation may be used for the evaluation of Z. Note that 

(3.5) the properties i), ii) also allow for a direct construction 
of the Fokker-Planck equation without making use of the 
Langevin equation. 

To calculate the nucleation rate, we have to construct a 
stationary nonequilibrium probability density. To main- 
tain a constant probability flux over the saddle point 



(3.3) to yield the Langevin equation 



dtMx, t) = -J2 M^J^LJiJJ{x, t) + C^{x, t) (3.6) 



where C^i is Gaussian white noise with {Q) =0 and 
2a 



{Q{x,t)Qix',t')) = S,,Six-x')6it-t') (3.7) 



where (. . .) denotes the average with respect to the Gaus- 
sian noise distribution exp{—f]A/ (4q;) / dtdx J2i &} ^nd 
f3 — l/kgT. The dynamics of the probability distribu- 
tion functional £i['0(a;)] = (Yii x ^('^^i-'^t^ ~ i'ii^))) with 
ipi{x,t) a solution of ( |3.6| ) is governecS by the Fokker- 
Planck equation 



dtg[ip{x)] 



dx 



SJ, 



5ipi{x) ' 



The probability current is given by 

1 5 



Hjipjix) 



(3A5i}j{x) 



(3.8) 



g[^{x% (3.9) 



In ( |3.9D we have exploited the antisymmetry of the off- 
diagonal part of M . Note that the current is only defined 
up to a divergenceless term. If we demand in addition 
that the equilibrium density has vanishin g cu rrent, the 
representation ( |3.9| ) is unique. Equations (3.S) and ( |3.9| ) 
have the following properties: 

i) The equations of motion for the thermal expecta- 
tion values (</5),(p) are identical to the expectation value 
of ( |3.3| ). [ Expectation values are defined by (1^) = 
J VipVp if g, where J "Dip denotes functional integration.] 

ii) The equilibrium density near the saddle point 



(3.10) 



with £s^^ as in ( 2.12D , is a stationary solution of ( |3.§| ) 
with vanishing current. Z is a normalization constant 



we impose the boundary conditions g ~ geq near the 
metastable state and ~ beyond the saddle point. 
Note that the realization of equilibrium at the metastable 
state requires a barrier energy which should be large com- 
pared to thermal energies. As a criterion we may use 
j3A£s^5. Since the prefactor is roughly of the order of 
the precession frequency 2jKe/MQ ~ 10^°s^^, this in- 
equality is satisfied even for very large switching rates 
and thus does not represent a restriction. The total rate 
is then obtained as the probability flux integrated across 
a surface transversal to the unstable mode. The deriva- 
tion is similar to that of LangeiQ and is presented in 
detail in the appendix. The switching probability per 
unit time of a particle with magnetization prepared in 
the metastable state into the stable state is then given 
by 



(3.11) 



where A is the cross sectional area of the sample and Eg is 
the energy (2.15) per area of the nucleus. The prefactor 
is given as follows 



PA detW'c /detH™P 



27r3 V det' m'-^l V detH 



(3.12) 



In ( 3.12| ) an explicit factor of 2 has been included, since 
the metastable state = tt, = 7r/2 may decay via ei- 
ther one of the two equi valen t saddle points (f)± . The first 
factor on the r.h.s. of ( 3.12 ) is the escape frequency of 
the unstable mode as defined by (3.2). This is the only 
term in (3.12) in which dynamical details of the system 
enter. The second factor arises from the integration over 
the zero frequency (Goldstone) mode and is defined by 
C — ^/£^L, where L is the system length and Sg is the 
energy per unit area (2.15). The third factor is also due 
to the Goldstone mode and determines the temperature 
dependence of the prefactor. The remaining factors ba- 
sically arise from the functional integration of the parti- 
tion function within Gaussian approximation ( 2.12] ) and 
(^.22| ). The determinants are defined as the products of 
eigenvalues 



detH™^ 



(3.13) 



where k denotes scattering states. The prime on det de- 
notes omission of the zero energy El'^ . However, note 
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that the unstable mode Eq"^ < enters ( ^.15 ) as if it 
were a stable one. The determinant of out of easy-plane 
fluctuations is defined as 



detH"P 



(3.14) 



The bound state energies Eq"^ and E^^ are given by ( 2.19| ) 
and (|2.20D. In paper I we obtained E^"^ ~ (5"^ (1,6.13) 



and Eq"^ and E^^ have been evaluated numerically for 
arbitrary R. 

Therefore we are left with the task of evaluating the 
products of the continuum eigenvalues. The continuum 
eigenvalues of Ti^'^ and TC^p coincide with those of 7^'""^ 
and 7i™P, respectively, and are given by 



Q 



2 

-1 



(3.15) 



Note that these equalities do not im ply a cancellation 
of numerators and denominators in ( 3.13| ) and (3.14), 
since the allowed k -values are different and fixed by the 
bondary conditions which we choose as periodic. 

The next sections are devot ed to th e explicit evaluation 
of the infinite products ( 3.13 ), ( 3.14 ) and the calculation 
of the escape frequency A+. 



IV. EVALUATION OF THE STATISTICAL 
PREFACTOR 



in the limits i? — > and R ^ oo differs from those for 
finite values of R. This casts some doubts on the useful- 
ness of suc h opera tors for an approximation of the prod- 
ucts ( 3.13 ), ( 3.14 ). However, by a careful analysis using 
the explicit form of Levinson's theorem in ID, we show 
that the exact fluctuation determinants converge to those 
evaluated by means of the limiting operators. 

The second method uses the explicit knowledge of the 
zero mode for the evaluation of the fluctuation determi- 
nants. This allows an exact determination of the statis- 
tical prefactor for the (/?- fluctuations. For p- fluctuations 
an analytical treatment is only possible in the limit 

These results can then be combined to obtain analyt- 
ical expressions for the total statistical prefactor in the 
limit of small and large nuclei as well as in the limits 



Q- 



and Q 



oo. In the limit Q 6 



U2 



corresponding to either small nuclei or large hard-axis 
anisotropy, the out of easy-plane fluctuations are sup- 
pressed and do not contribute at all. While the latter 
result is to be expected from the fact that out of easy- 
plane fluctuations are suppressed due to their mass Q~^, 
the former result is somewhat surprising. It is related 
to the divergence of the characteristic length scale for 
R ^ which renders even a small hard-axis anisotropy 
effectively large. In both limits the system may be de- 
scribed by an effective model discarding the out of easy- 
plane degree of freedom. As we will later address, this 
model is equivalent to a double sine-Gordon model in the 
azimuthal variable 6. 



In the following we describe two meth od s for the eval- 
uation of the statistical prefactors ( 3.13 ), ( |3.14 ). 

The first method (see e.g. Ref. p5[) is based on the 
knowledge of scattering phase shifts of the continuum 
eigenf unctions as well as the bound state energies of the 
operators H'^'^, H^p. In strictly one dimensional prob- 
lems one has to distinguish between the scattering phase 
shifts of odd and even parity wavefunctions. This is in 
sharp contrast to the familiar situation of a 3D s-wave 
scattering problem where the scaled wavefunction indeed 
obeys a ID Schrodinger equation but is required to van- 
ish at the origin. Surprisingly, this issue has been ignored 
until rcccntlyliilllj, possibly also due to the fact that the 
most widely used ID potentials belong to the rather spe- 
cial class of reflectionless potentials for which scattering 
phase shifts of even and odd wavefunctions coincide. 

So far, the present method for the, evaluation of the 
prefactor has indeed been used onlytj in the case of re- 
flectionless potentials where the scattering phase shifts 
of even and odd wavefunctions coincide. Here, however, 
the potentials in Ti*''' and T-C^p only become reflectionless 
in the limits R —^ and i? — )■ oo. For intermediate values 
of R, the corresponding potentials are not reflectionless 
and scattering phase shifts of even and odd parity wave- 
functions have to be distinguished. 

There is another surprising feature of these operators. 
The number of bound states of the operators which arise 



A. Scattering phase shift method 



In this section we evaluate the products (3.13), (3.14) 
using the knowledge of bound state energies and scatter- 
ing phase shifts of the operators H'^'^, which have 
been evaluated in paper I. 

In order to evaluate the density of states we con- 
sider eigenfunctions obeying periodic boundary condi- 
tions. The modes around the metastable state (0m, ^m) — 
(tt, 7r/2), are then the pla ne w ave eigenfunction eigen- 
functions of the operators ( 2.23 ) which can also be writ- 
ten as sinkx, coskx with wavenumbers 



k = 



2nn 



(4.1) 



where n = 0, 1, . . . for even parity and n = 1, 2, . . . for odd 
parity continuum eigenstates. The corresponding density 
of states is 



mi 



dn 
Ik 



L 

2^' 



(4.2) 



At the saddle point {(f)s,ds) = (0s(a;), 7r/2), however, 
we encounter a different situation. The nonuniformity of 
the nucleus, i.e. the nonconstant potentials in 7i*'^,?i'*^' 
lead to phase shifts of the continuum eigenfunctions. In 
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contrast to 3D problems, where the wavefunction always 
vanishes at the origin, we have to distinguish between the 
phase shifts of even (e) and odd (o) parity wavefunctions. 
We define phase shifts as in paper I 

Xl%) {x ^ ±00) (X cos (fcx ± (fc)/2) , (4.3) 
Xl^^^ix ^ ±(X3) « sin (fcx ± A|„)(fc)/2) , t^^,p. (4.4) 

Since only the eigenfunctions of H." "^'^ exhibit a phase 
shift, we have omitted the superscript s on A. Note that 
all phase shifts also depend on the p aram et er R . Periodic 
boundary conditions together with ( |4.3| ), (4.4) imply 



fcL + A|^-)(fc) = 27rn, 



1,2,.. 



(4.5) 



where i — ip,p and j — e,o. EoUowing the arguments 
of Ref. |Tl|, the lowest allowedEj fc-values in (4.5) are 
k ~ lixjL for odd parity eigenfunctions and k — ixjL 
for even parity eigenfunctions. Note the surprising fact 
that the latter value does not coincide with the lowest 
fc- value ( |4.1| ) of even-parity solutions in the absence of a 
potential. 

The density of states for odd-parity continuum eigen- 
f mictions follows from (4.5) 



dn 
dk 



L 
2^ 



1 dAl^(k) 
27r dk 



(4.6) 



Since the spectrum of even-parity continuum eigenfunc- 
tions starts at k — tt/L while free solutions start at fc = 0, 
the density of states exhibits an additional J-function 
contribution at A: — s- 



L 1 dAUk) 





27r dk 



Isik^ 



0+), 



(p,p. 



(4.7) 



We are no w in a position to express the ratio of the 
products in ( 3.13 ) in terms of the density of states: 



= cxp 



00 

dk 



xlni^r}, (4.11) 



where i — (p,p. After inserting (|3.14| ), ( [4. 61) , (|4.7| ) into 

(4.11), using (4.1c) and performing a partial integration 
we obtain for the 95-fluctuations 



U,' Elf 



exp 



dk r 
2^ 



Af, 

(e) 



A^ 

(o) 



2kS^ 



k^S^ 
(4.12 



where we used the fact that the phase shifts vanish as 
l/k6 for fc — > cxD according to Born's approximation 
(1,6.29). In a completely analogous way we obtain for 
the p-fluctuations 



U,'El^' 



X exp / — 



dk r 

2^ 



(e) 



2k6^ 



1 + 0-1(52 + fc2J2 



(4.13) 



Note that in contrast to (4.11), the integrands in ( 4.12| ) 
and ( 4.13 ) are independent of the fc — s- limit which due 
to Levinson's theorem is se nsitiv e to the num ber of bound 
states. This fact renders ( 4.12 ) and ( 4.13| ) suitable for 
phase shift approximations that converge nonuniformly 
to the exact phase shifts for fc ^ 0. 

In the next two subsections we explicitly evaluate the 
prefactor in the limit of small and large nuclei. We show 
that taking the limit i? — > 0,oo of (4.12). ( 4.13D is equiv- 
alent to a direct evaluation of the determinants of the 
operators that arise in these limits. 



This 5-function also ensures that the number of states of 
the free problem equals that of the scattering problem 
including bound states 



dk 



0)' 



(4.8) 



1. Prefactor for h ^ 1 

In the limit i? — > 0, when the external field is very close 
to the anisotropy field, the nucleus represents a slight but 
spatially extended deviation from the metastable state. 
The operators H^'^ and H'^ given by ( |2.13D and ( |2.14| ) 
then reduce to 



where, according to (1,6.25) 



1, 



(4.9) 



are the number of even - and odd-parity bound stat es o f 
and W^. Eq (|!|) wfth (Q is verified using (0), 
(4.7) and ID Levinson's theorem (1.6.23,24) which states 
that 



Af,)(fc = 0)=^, AP^)(fc = 0) = 2^, 
Af^^(fc = 0) = 3^, A^^^(fc = 0) = 27r. 



(4.10) 



7Y^^ = -^+r2(l-6sech2^) 



da; 2 



+ (5^2(1 - 2sech2 



(4.14) 
(4.15) 



The potentials in ( 4.14 ), ( 4.15 ) are refiectionless and the 
solution of the corresponding eigenvalue problems is well 
known (see also the appendix of paper I). There are 
bound states with the following energies 



El"^ = 0, 



(4.16) 
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Note that the two eigenvalues El\ of and 
turn into zero energy resonance£3 of Ti,^ '^ and H,^^ and 
have therefore no c ount erpart s in ( 4.1(j ). The contin- 



uum eigenvalues of ( 4.14 ) and ( 4.15 ) are given by ( 3.15 ) , 
respectively. Since the potentials are reflectionless the 
scattering phase shifts are parity independent and given 
by 



A'^(fc) = 2arctan 



(M)2 -2' 



AP(fc) = 2arctan 



1 



and their long wavelength behavior 

A'^(fc 0) = 27r, AP(fc^0)=7r, 



(4.17) 



(4.18) 



(4.19) 



detH""^ 
de7|H^ 

det mp 



(4.24) 



For the evaluation of the r.h.s . in ( 4.24|) we proceed simi- 
larly to the derivation of ( 4.12 ) and ( [4.13 ) with the follow- 
ing modifications: The density of states p^'t^ — p^*^ (fc) in 

( PH) has to be replaced by Pfj)i^) = -^dA'/dk 

with A* given by ( [4.17| ) and ( |4.18D {i ^ (p,p). Usi ng the 
versi on ( 4A9) of Levinson's theore m tog ether with (4.17) 
and ( 4.18 ), we recover the results ( 4.22 ) and ( 4.23 ) after 
integration. 

To summarize, we thus have shown that the small nu- 
cleus approximation may be used for the evaluation of 
the fluctuation determinants, or explicitly 



is in accordance with Levinson's theorem (1,6.27) for re- 
flectionless potentials. As has been discussed in paper I, 
Sec. VLB, the convergence of A''', A'' towards the ex- 
act phase shifts A^*^ is in general nonuniform for fc = 

(cL Figs. L6,7). However, since the integrand in (4.12) 
and ( 4.13 ) vanis hes fo r k = 0, w e can safel y inser t the 



approximations ( |4.17D and ( |4.1^ ) into ( |4.12| ), ( |4.13| ), re- 
spectively, and obtain 



lim 



UkEr 



Uk' El 



36 r 



and 



lim 



u,Er 



Elf 



(4.20) 



(4.21) 



respectively. For the evaluat ion o f the determi nant s, we 
substitute the bound states (4.16) for those in (3.13) and 
( |3.14 ). However, we have to complete (4.16) by the zero 
energy reson ance s £2"^ — S~ ^, E^^ = Q^^ + 6~ 
gether with ( [4.20D and ( [4.21D , we then obtain 



To- 



and 



hm r- = 12 6 ^, 



hm — TT, — ~ li 



(4.22) 



(4.23) 



where in (|420| ), (g^, 6-'^ = . The resuh (|423 ) 
is remarkable since it shows that the fluctuations in p- 
direction do not play a role at all for small nuclei, inde- 
pendent of the size of the hard-axis anisotropy constant. 
This suggests that in the limit i? ^ the system may ef- 
fectively be described by a double sine-Gordon equation 
in the azimuthal angle 4>. We shall return to this issue in 
Sec. VH. 

Alternatively, although less carefully, we can inter- 
change the limit i? ^ with the functional integration, 
and directly calculate 



detH™^ detrr-'P 
lim — — 

B.^Q det' \H''v\ det'lH^-^l' 

det7^'"P det7^'"P 
hm — = = — . 

i?-^o deiWP deiWP 



(4.25) 



where the r.h.s. has b een evaluated in leading order in 
R. The relation (4.25) has not been clear on the onset, 
since the operators Ti.'^'f ^ T-C^p exhibit a different number 
of bound states and different long wavelength behavior of 
the scattering phase shifts than the operators Ti^"^ , Ti^P . 
As we have shown now, these two differences co nspire in 
such a way that an interchange of the limits in (4.25) is 
indeed correct. 



2. Prefactor for h 

For i? ^ 1 the nucleus separates into two independent 
TT-Bloch walls. Correspondingly Ti*''' and merge into 
the same operator Ti" which consists of two independent 
potential wells of the form —2S~'^sech^{x/6 ± R). The 
bound states of Ti.'' are then given by the symmetric and 
antisymmetric linear combinations of the ground states 
of the single wells and have energies (cL (1,6.4,6.6)) 



rpsp 



'2R 



E'/ = 0, 



8e 



-2R 



(4.26) 



Note, that in this approximation £'2''' has merged into a 
zero energy resonance of 7i*. The continuum eigenvalues 
are identical to (??) while the scattering phase shifts are 
twice those of a single potential well (cf. ([4.18 )) 



A" (A;) = 4arctan— . 

ko 



(4.27) 



The coincidence of the phase shifts of even and odd eigen- 
functions originates from the fact that the two —2 sech^ x 



potential wells are reflectionless. The phase shifts (4.27) 
obey the reflectionless version of Levinson's theorem, i.e. 
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Delta'' {k ^ 0) = ttN with N the number of bound 
states. But as in the previous subsection, the phase 
shifts A'^j and A^^^ only converge on the open inter- 
val < fc < oo towards A". Inserting (4.27) into 

^= ,5-2 



( 3.14 ), we obtain 



( 14.121), (^.13|), a nd using ( f4.26| ) together with E: 
in (13.1^) 



detH™^ ^2e2« 



lim 

ii^oo det' iH'^fi 



(4.28) 



slice representation of the path integrals. This method 
allows for the exact evaluation of the statistical prefactor 
in the yj- variable for all values of R. In the limit — 
we are also able to evaluate the prefactor for out of easy- 
plane fluctuations. We first show how this method can 



be applied to the evaluation of ( p.l3| ) . According to Ref . 
p3 we have 



deti 



D^°\l/2) 



(4.32) 



lim . iVT 



Q + VQ 



(4.29) 



The latter result approaches 1 for large hard-axis 
anisotropies as expected. However, in the opposite limi t 
of high Q, the p-fluctuations lead to a prefactor (3.12) 
proportional to Q. 



3. Out of easy plane fluctuations for Q — > oo 

A large hard-axis anisotropy leads to the suppression 
of out of easy-plane fluctuations by the existence of a 
large mass gap. Therefore we expect the fluctuation de- 
terminant to become one in this limit. 

(o) 



To prove this conjecture, we remark that A^^^. 



are continuous functions which are proportional to 1/kS 
for fc — > oo according to Born's approximation (1,6.29) 
and remain finite for fc due to Levinson's theorem. 
Therefore both phase shifts obey the inequality A(fc) < 
c/k S wit h a suitably chosen constant c. For the integral 
in ( |4.13| ) we thus obtain the following inequality 



0< — 



dk 

TT 



kS^ 



(o) 



< 



(4.30) 



where d = (c^^^ -|-c^^j)/2. Since d is independent of Q, the 
upper limit tends to zero for Q^^S^ — s- 00 and therefore 



we have with E^^ = Q 



e < 1 and (4.13), (3.14) 



E'/ 



Q-i + eS~^ with < 



lim 



detH"P 



= 1. 



In fact we have now shed some new light on the re- 
sult (4.23). The hard-axis anisotropy does not enter in 
isolated form but rather in the combination Q~^S^ — 
coth^ R which shows that due to the diverging length 
scale, the hard-axis anisotropy becomes effectively strong 
for small R, no matter how small is. 



B. Jacobi method 

There is also an alternative way for evaluating the 
products of eigenvalues which has its origin in the space 



The notation det^ on the l.h.s. of ( 4.32| ) indicates that 
the evaluation of the determinants relies on eigenvalue 
problems defined on the finite interval [— L/2,L/2] with 
respect to functions that vanish at the end of the interval. 



The functions and on the r.h.s. of (4.32) are 



yip ciiiu. j^ip 

solutions of the differential equations 



dx'^ 



n'^^D^ix) = 0, 
+ 5-Ad(^\x)^0, 



(4.33) 
(4.34) 



with the following "initial" conditions (the prime denotes 
d/dx) 

D^{-L/2) = 0, D'^{-L/2) = \ (4.35) 
D^°\-L/2) = Q, Df{-L/2)^l. (4.36) 

Note that on the finite interval the first excited eigen- 
function of Ti'"^ has no longer zero energy and therefore 
the l.h.s of (4.32) is well defined. The eigenvalue problem 
of this quasi zero-energy mode is written as 



(4.37) 



where for large system lengths L, /i > is small, and / 
obeys the boundary conditions /(±i/2) —f'{±L/2) = 0. 
Note that for L ^ 00 we have / Xi^ and fj. ^ 0. The 
fiuctuation determinant is then obtained as follows: 



det' m^^l 



lim 

L — 'OO 



dI'\l/2) 



D^{L/2) 



(4.38) 



(4.31) We now turn to the evaluation of the r.h.s. of (4.38) 



D^\l/2) is easily obtained by integration of (|4.34|) with 



(4.36) 



D^^^ (L/2) = S sinHL/S) ^ -e^/^ (4.39) 

Th e fun ction D^{x) obeys the sa me d ifferential equa- 
tion ( 4.35 ) as the zero mode x^' ^(x) ( 2.20 ), but subject to 



different boundary conditions (4.35). Therefore Dip(x) is 
a linear combination of x'l'ix) and the unknown linear 
independent solution ^^'''(x) of the differential equation 
(p3|) 



D^{x) 



<r(^), 



(4.40) 
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with real constants u, v. Xi^ is taken to be unnormalized 



-P = [sech(a;/^ ~ R) - sech{x/S + R)] 



dx 



(4.41) 



The normahzation of ^l'^{x) is chosen such that the 
Wronskian 



Sip 

Xi 



dx 



dxZ 
1 dx 



Sip 



= 1. 



(4.42) 



In order to satisfy the initial conditions (4.35), we must 
have 



-er(-i/2), 



xri-L/2). 



(4.43) 



From (4.41) we infer the asymptotic behaviour 

^ iVsgn(2;)e"l"^l/^ for x ^ ±oo, (4.44) 



with 



N = 4(5^1 sinhi?. 



(4.45) 



The symmetry of the potential in Ti.'^'^ and the antisym- 
metry of Xi'^ allow us to choose ^J"^ as a symmetric func- 
tion which has the asymptotic behaviour 



e^"^ -> iV'el^l/^ for X -> ±oo. 



From (|4.42D it follows that 



N' = — , 



and therefore with ( [4.43D -( ^!47| ) we have 



DAL/2) 



(4.46) 



(4.47) 



(4.48) 



Finally we are left with the evaluation of fi. The eigen- 
function / may to first order in fi be expressed as 



L/2 

f{x) = T]{x) + ^ / G{x, y)r]{y)dy, 

-L/2 



(4.49) 



with the Green's function 

G{x,y) = e{x - y)[xr{x)^r{y) - xrivKn^)]- (4.50) 

The quasi zero-energy eigenvalue fi is now determined 
such that /(±L/2) = 0. The func tion 77 is a solution 
of the homogeneous problem ( 4.33 ) which satisfies the 
boundary condition 7]{~L/2) = 0, 



with 



c=-xr(-L/2)/er(-v2). 

The requirement /(i/2) = then leads to 



(4.51) 
(4.52) 



A* = 



x7{L/2)+cir{Ll2) 



j%dy[{xr{y)nr{L/2) 



cxriL/2){cr{ym 

(4.53) 

Since the normalization of x'l' is independent of L, the 
first term in the denominator is of the order exp(L/2S) 
whereas the second vanishes as exp(— L/2(5) and thus may 
be neglected. In leading order in exp{~L/S) we thus 
obtain 



2xiiL/2) 



UL/2)J%xUy) 



64 (5-3 ginh^ R 
£JR) ' 



-L/S 



(4.54) 



where w e hav e mad e use of ( 4.41|) a nd (1,3.12). Inserting 
( fl.39| ), ( [4.48[ ) and ( [4.54[ ) into ( [4. 38] ) and performing the 
limit L 00, we finally obtain the result 



detTi""^ 8tanh3i?sinh^i? 



""^ tanhi?- i?sech^i? 



det' \n 



(4.55) 



Note 



which is exact for all values of the extern al fie ld, 
that in the lim its R ^ an d R 00, Eq. (4.55) reduces 
to the results ( [4.22|) , (|4.28| ), respectively. 



The above metho d can not be used for the evaluation of 
the p-determinants ( 4.58| ) for arbitrary values of the hard- 
axis anisotropy since the zero energy eigenfunction of Ti"'' 
is not explicitly known. In the limit of small Q~^S^, 
however, the fluctuation determinant may be calculated 
exactly: 

The Q-^-independent operator H^p ~Q~^ exhibits the 
zerOjenergy mode Xa' we can proceed along the same 
linegl3 as in the derivation of ( 4.55| ) to obtain 



8 tanh^ R cosh^ R 



det(H"P - Q-^) 

det'iWP - g-i) ~ tanhi?-hi?sech^i?' 



(4.56) 



where the prime d enote s omiss ion of the zero en- 
ergy mode. From ( |4.11 ) with ( 3.15[ ) it follows that 

0{Q ^(5^). With ( ^.14[ ) we then obtain in leading order 
in the small parameter Q^^S^ 



detH^P 
det WP 



1 



Q-i det(H™P - Q- 



^0 



E 



det'in^P - Q- 



E'P -Q-^ 8 tanh^ R cosh^ R 
Ef tanhi? + i?sech^i?' 



(4.57) 



The second factor on the r.h.s has been retained since 
the coefficients of its power expansion in Q~^5^ diverge 
for R —> CXI. However, in the limit Q^^ ^ E'^ (R not 
too large), it reduces to one. Note that the result (4.57) 
is only valid for Q~^6^ small and it can therefore not be 



used for R small. In this latte r case ( 4.31 ) ap plies. 

With the exact result ( 4.55| ), the prefactor ( [3.12] ) now 
takes the following form 



4 / det T-C^P 

n = A+LV/Jyl-^tanh^/^i? sinhj?^ , (4-58) 



which is further evaluated in the following section 
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V. NUCLEATION RATES IN THE 
OVERDAMPED LIMIT 

In this s ectio n we shall derive analytical results for the 
prefactor (4.58) in the limit of large and small values of 
Q~^S^ as well as large R. In the intermediate parameter 
range the prefactor is evaluated numerically. The discus- 
sion in this section is restricted to the regime of large 
damping . The case of moderate damping which is more 
relevant for real systems shall be discussed in the next 
section. 

We start with the evaluation of the decay frequency A+ 
of the nucleus. Fo r la rge values of the damping constant 
a, the equations (3.2) characterizing the unstable mode 
of the nucleus decouple and take the form 



X+P+ = -aWPp+, 



(5.1) 



with A+ > 0. The unstable mode is thus given by the 
ground state of , i.e. ((/3+,p+) cx (xo'^'jO) ^"^^ for 
large a the corresponding escape frequency is given by 



A. 



{R)\. 



(5.2) 



The eigenvalue Eq^{R) has been investigated in paper 
I . It is know n analytically in the limits i? — > 0,cx) (cf. 
(4.16), ( 4.26| ) respectively) which allows an explicit eval- 
uation of the prefactor in these limits. 

i) For large nuclei (i? oo) we may use ( 4.26|), (4.29) 
and together with £s 
the following form 



4 and the prefactor ( 4.58 ) takes 



16 

73/2 



(5.3) 



u) For Q~^(5^ > 1, ([4.3l[) and ( |5^) may be inserted 
into (|4.58D to yield 



aLJl3A^rj,T\E''^{R)\t8inh.^/^R sinhi?. (5.4) 



For small R, this reduces to 



n = aL./l3A^R^'^. 



(5.5) 



Note that the results, (x4) and ( p.5| ) are both indepen- 
dent of the value of the hard-axis anisot ropy. 



ui) In the limit Q^^S^ ^ we can use ( |4.58|) and ( |4.57| ) 
to obtain 



r3/2 I 




E{^ - Q-i 
E{^ 



tanh^ R sinh^ R 
\/ tanh i? + i?scch^i? 



(5.6) 



For values of R such that Q ^ <C i?*^ (i?) , the square root 
cont aining E''^ r edu ces to 1. In the limit R oo, Eqn 
(Eel) reduces to (O) with Q-^ ^ 



Note however, that the present theory is only valid 
for hard-axis anisotropics which are not too small such 
that the amplitude of out of easy-plane fluctuations is 
much smaller than one. Requiring the thermal expec- 
tation value (p^) with respect to the Boltzmann weight 



e^p{~l3A£r'} with Sr' as in ( |2.12| ) to be smaller than 
one and noting that the lowest eigenvalue in p-direction 
is Eq^ — Q~^, we obtain for the validity of the present 
theory the condition (3A\fAKlKh/ K^. > 1. 

For all other parameter values the prefactor has been 
evaluated numerically. The corresponding results are 
shown in Fig. | with the dashed and d ashe d-dott ed l ines 
representing the asymptotic form ulas ( ^.3[ ) and ( |5.5| ) re- 
spectively. Re insta ting the units ( ^.4|) we recognize that 
the prefactor ( 3.12| ) is inversely proportional to p, a fact 
which is in accordance with the general behaviorllj of nu- 
cleation rates in the overdamped limit. 



VI. NUCLEATION RATE FOR MODERATE 
DAMPING 

In the last section we have presented results for nu- 
cleation rates of thermally activated magnetization re- 
versal in the overdamped limit. In that case, the decay 
of the nucleus is governed by purely dissipative mecha- 
nisms. However, in ferromagnetic materials the damping 
constant is usually of the order a ~ 10~^ or sometimes 
as small as 10^"* in some high purity materials such as 
YIG. 



According to (3.12), the dynamicr^roperties of the sys- 
tem enter the nucleation rate onlyH in the form of the 
decay frequency of the nucleus. In order to evaluate 
the nucleation rate in the moderately damped regime we 
therefore have to include the conservative, precessional 
part of the dynamics for the evaluation of the decay fre- 
quency A+. In the limits of small and large nuclei this 
decay frequency can again be expressed in closed analyti- 
cal form thus enabling us to give exact results of the total 
prefactor and hence of the total rate. To the best of my 
knowledge, this, provides the first application of Langer's 
general theoryQ of moderate friction to a system with 
infinitely many degree of freedoms. 

In a first part we discuss the escape frequency A+. 
We obtain exact expressions in the limits R ^ oo and 
Q~^d^ —> oo as well as an approximate formula which 
expresses A+ by Eq"^ . These results allow an exact eval- 
uation of the nucleation rate in the limits R and 
R — > oo. 



A. Escape Frequency 

The dynamically unstable mode {(p+,p+) of the nu- 
cleus is the solutio n of the (non-Hermitian) coupled 
eigenvalue problem ( |3.2| ). Before turning to a quanti- 
tative analysis we give a qualitative discussion of the 
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parity and relative sign of the functions (/9+, p+. Both 
functions are nodeless and symmetric in x with opposite 
relative sign as may be seen from the followi ng p lausi- 



bility arguments: For a ^ oo we know from (5.1) that 



the dynamical unstable mode coincides with the ground 
state of H'^'-f, i.e. (93+, p+) cx (xo''',0), and therefore (p+ 
is a symmetric nodeless function in x. For finite values 
of a there will be a nonzero p_|_-component. Since the 
nucleus represents an untwisted 7r-Bloch wall pair (cf. (I, 
3.11 )), the instability represents a confluence or a sepa- 
ration of the two domain walls which is associated with 
a monotonical increase or decrease of the angle 0. Hence 
ip^ is a symmetric nodeless function in x. To comment 
on we have to recall that a motion of the domain 
wall is only possible it-the structure exhibits an out of 
easy-plane componentc3. In order for the domain walls 
to move in opposite directions, the out of easy-plane com- 
ponent must have the same sign at the center of the two 
oppositely twisted kinks and due to the gyroscopic na- 
ture of the equations of motion, p+ and must have 
opposite signs. 

Therefor e we are looking for even-parity nodeless so- 
lutions of (3.2) with opposite signs. The ambiguity in 
the overall sign of ((^+,p+) describes the freedom of the 
nucleus either to collapse or to expand. Inspecting (3.2) 
we recognize that the eigenvalue problem can be easily 
solved if Lp+ and p+ are the ground states of Ti.'^'^ and H'^p 
respectively and proportional to each other. This is ful- 
filled in two limiting cases, i) R large a nd ii) Q~^S^ ^ 1. 

i) For large R we have according to ( ^.19 ) and (I, 6.2) 



Xo oc Xo ^ sech(x/(5 + R) + sech{x/S — R). (6.1) 
Inserting oc p+ cx Xa^ into (3.2) and using that the 



ground state energy of H.'^'p is given by 
obtain 



A- 



/(I) 



WG 



(6.2) 



where E^ 



-8e~^^ in the limit i? — > 00. The square 
root in (6.2) has been retained because the relative mag- 
nitude of the (small) parameters E^'^ , a, and has not 
been specified yet. The corresponding unstable mode is 
given by 



^,p+) oc |^sech(- 



R) 

X (1,-Q[A+ 



sech(- 




R) 



aE'/]). 



(6.3) 



The plus sign of the square rooot is chose n in order to 
reproduce the correct asymptotic behavior (5.2) for large 
a. Note that Eq. (|6.3| ) agrees with the statements made 
above: The functions are symmetric and nodeless 

while the ratio p+ / 1^9+ is always negative and vanishes for 
a ^ 00. For a ^ we have p+/ip+ — — ^/Q\E^^{R)\ . 

ii) For Q-M^ > 1 we have n'Pp+ = Q-^p+ + 0{S-^). 
Thus the first Eq. of (3.2) can be solved for p+ and 



after insertion into the second eq. of (|3 
eigenvalue problem is obtained: 



= -X, 



X+Q 



1 + + aX+Q 



the following 



(6.4) 



The solution of this equation is known, i.e. ip^ cx Xo''': 
and hence the coefficient of (p+ on the r.h.s. equ als E^''^ . 
Solvi ng for A+ we recover the expressions ( |6.2D for A^ 



and (6.3) for p^/ip^ but with Eq^{R) now evaluated for 



arbitrary values of R. This i s a remarkable result as it 
demonstrates the validity of (6^), (|6 
limits Q-^S^ > 1, i? > 1. 



in the opposite 
Note that for large Q^^, 
and (6^) hold for all values of R. In the particu- 
lar case of small R {6^ ^ 1) we can inser t th e small R 
approximation E'^'f = -SR'^ + 0{R'^) into (lej) to obtain 



A, 



- 3i?2 



[Q-2 +60-ii?2] _^3g-i^2^ (6 5) 



The unstable mode is then given by 

iip+,P+) cx sech2(5)(l, -g(A+ - 3aR^)). 



(6.6) 



The validity of the expression ( |6.2| ) in the opposite limits 
i? — > 0,cx) might hint to a more extended validity. In 
order to investigate A+ for intermediate values of R at 
arbitrary we have to resort to numerical methods. 



It turns out that the direct integration of (3.2) yields 
rather inaccurate results (errors of 10% ). Considerable 



improvement has been achieved by converting (3.2) into 
two decoupled 4**^ order differential equations in ea ch o f 
the variables ip+,p+. As is seen from Fig. ^, Eq. (6.2) 
provides an excellent approximation to these numerical 
results. 



B. Nucleation Rates 

We are now in a position to give the results for the 
prefactor for moderate damping. Results in closed form 
are obtained in the limits R ^ and i? — > 00. For 
Q~^S^ —^ 00 the prefactor can be expressed in terms of 
the negative eigenvalue Eq^ which for arbitrary values of 
R has to be evaluated numerically. 



i) For large R, we can combine ( [l.58| ), ([4.29[) and (|6.2| ) 
with E^"^ = -8e-2^ to obtain 



(6.7) 




[Q-2 + 16Q-ie 



The square root has been retained since the relative order 
of the small parameters a and e~^^ has not been specified 
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yet. However, in an expansion of the square root only 
leading terms in e~'^^ should be taken into a ccoun t, 
ii) For Q-i(52 _^ oo, the results and yield 



= tanh^/2 R sinhi? 



(6.8) 



with A+ given by (S^). For small R, this result reduces 
with (^.5|) to 



in the limit Q~^5^ oo (i.e., large hard-axis anisotropy 
and/or fields close to the anisotropy field), magnetization 
reversal rates become equivalent to the creation rates of 
kink-antikink pairs in the double sine-Gordon model. 

Discarding noise terms for the moment, we consider 
the dynamics of a field variable (t>{x,t) which is governed 
by the damped double sine-Gordon equation 



Qd^(j) + adt(j) = - 



S£, 



dSG 



(7.1) 



n = lJ(3A 



T-3/2 



R 



.5/2 



with the energy 



/ax 2 1 fdSG = / dx < -(dx(j))^ + -sin^ (h- hcoa(j)> . (7.2) 

(f) [Q-' + 6Q-ii?2]+3Q-ii?2l. (6.9) Jl/2 [2' ''^' 2 ^ -A-J V ; 



Also here, the square root has been retained since we did 
not specify the relative magnitude of the small parame- 
ters a and R. However, we have to keep in mind that 
upon expansion of the square root only terms in leading 
order in R have to be kept in order to be consistent with 
the derivation of the determinant of p- fluctuations. For 
small damping constants a ^ y/QR, Eq. ( |6.9| ) reduces 
to 



4V3 



n = l^pa^Vq^r^/^ 



(6.10) 



This limit is realized in typical experimental situations. 
For y/QR <C a we obtain 
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a . 



(6.11) 



For large values of a, Eq. (p. 11) merges into the prefactor 
for the overdamped regime ( |5.5| ). 

The above results for the prefactor are summarized in 
table I. In Fig. 5, numerical results for the prefactor are 
shown for arbitrary values of R. The prefactor is maximal 
for i? ~ 1 and decreases as the external fleld approaches 
the anisotropy field, i.e h ^ 1, or as the field approaches 
zero. One clearly recognizes that fl is independent of the 
hard-axis anisotropy and the damping constant a 
for small and intermediate values of R, respectively as 
predicted by (|6T|), (|6T0| ). 

T he to tal rate for magnetization reversal is then given 
by (3.11) and experimental consequences of this result 
shall be discussed in section VHI. A detailed discussion 
of experimental implications of these results may also be 
found in Ref. |42, 



VII. 



RELATION TO THE DOUBLE 
SINE-GORDON SYSTEM 



In this section it is shown that the results of the pre- 
vious sections allow us to calculate the nucleation rate 
of kink-antikink pairs of the double sine-Gordon model 
for moderate to large friction. Second, we shall see that 



The constant Q plays the role of a mass and a is a 
damping const ant . In the overdamped limit, the inertia 
term Qd^cj) in (7J) can be neglected and the dynami cs i s 
purely determined by the damping term. Note that ( |7.2| ) 
is equivalent to the energy density (|2.5|) restricted to the 



easy-plane 9 = tt/2. Therefor e, E q. ([7.l[) exhibits the 



same saddle point solution (j)^ (2.9) as the full magnetic 
system. The corresponding barrier energy between the 
metastable state </> = tt and the absolute minimum cj) = 
is given by £s ( 2.15| ). Expanding 0(.t, t) ==^{x) + (p{x,t) 
and linearizing the equation of motion ( [7.l[ ) around the 
saddle point yields 

Qd^ip + adfip + H'''^ip = 0, (7.3) 

where W**^ is given by ( 2.13| ). A description of the 
stochastic dynamics in the vicinity of the saddle is ob- 
tained by adding the stochastic force ^{x,t) to the r.h.s. 
of (7.3). with the noise correlation {^{x' ,t')£^{x,t)) = 
{2a/P)6{x — x')S{t — t'). The corr esp onding Fokker- 
Planck equation takes the form of (3_^) if we identify 
f3 = [3A, Mil = 0, Mi2 = -A'f2i = 1, M22 = a and 

The transition rate from the metastable state (neglect- 
ing transitions that lead back from the metastable state 
over the barrier to the initial state) can then be calcu- 
lated as in the appendix with the result: 



dSG 




(7.4) 



where £s is given by dlTsI ). In (|7j) a factor of 2 has 
been included due to the existence of two equivalent sad- 
dle points ±(/)s, and the ratio of the determinants has 
been calculated in ( 4.55| ). A+ > is the nucleus decay 



frequency which is obtained by insertion of Lp 



into (7.3), 



l + 4Q|£;o^n/a2 



(7.5) 



with E]^ th e (n egative) ground state energy of TC^f . 



Eqns (7^), (7^) constitute the creation rate of kink- 
antikink pairs in the moderately damped double sine- 
Gordon system. 
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It may now be veri fied tha t the magnetization rever- 
sal rate ( 3.1l| ), ( |3.12 ) with ( 1.31 ) is equivalent to the 
result (7.4) in the Hmit Q~^5^ oo, provided that the 
time scale in (|!|)-(U) is chosen as [t] = Mo/{2'-fKe) 
while energies and lengths are chosen as in_(2.4). Taking 



the limit QlE^'^l/a^ in (|6^ and ([7j)^e obtain 
A+ = (a + a-i)|So'^| and A+ = l-Eg'^l/a, respectively. 
Reinstating units, the equivalence of F and T^sg is im- 
mediately verified. 



VIII. DISCUSSION 

In the previous sections we have investigated the rate 
of magnetization reversal in an effectively ID ferromag- 
net which describes magnetization configurations in an 
ideal elongated particle of a small constant cross section 
A. The experimentally most important conclusion is the 
existence of a saddle point structure which is localized 
along the sample. Unlike the Neel-Brown theoryo which 
leads to a barrier energy VKf,{l — h)^ proportional to 
the particle volume V, the present theory leads to an 
energy barrier A£s that is proportional to the sample 
cross section and t o th e domain wall energy (after rein- 
stating the units (2^)). For sufficiently elongated par- 
ticles the energy of the nonuniform barrier is thus al- 
ways lower than that of the uniform one and thus the 
present theory predicts much lower coercivities than the 
Neel-Brown theory. To illustrate this, we consider the 
following typical material parameters of particles such as 
Cr02: A = 5 ■ 10"'^ erg/cm, Ke = 7 • 10^ erg/cm^. 
Mo = 480 Oe, 7 = 1.5 • lO^Oe-^-'^. For T = 300 K, 
= 0. 2, a = 0.05, the numerically evaluated switch- 
ing rate (3.11), (4.58) is shown in Fig. |^ as a function 
of the external field for various particle diameters but 
for a fixed aspect ratio of 15:1. The dotted lines repre- 
sent the predictionia of the Neel-Brown theory, while 
HcxtMo/{2Ke) = 1 is the Stoner-Wohlfahrth value of 
the nuclcation field. One clearly recognizes the dramatic 
coercivity reduction for particles with small diameters. 
Note that the switching rate at a given field predicted by 
the present theory exceeds that of the Neel-Brown theory 
by more than 10 orders of magnitude. Conversely, the co- 
ercivity of a particle of diameter 100 A exhibits a coerciv- 
ity which (depending on the measurement time) is about 
one third of the Stoner-Wohlfahrth value. Since the bar- 
rier energy A£s is independent of the particle length, the 
present theory predicts the coercivity to become indepen- 
dent of the particle length for sufficiently long particles. 
This is in contrast to the theory of Neel and Brown which 
predicts a suppression of thermal effects in the particle 
volume. 

Experiments investigating the coercivity of a single 
elongated particle indeed show a significant coercivity 
reduction from the Stoner-Wohlfahrth value for fields 
along the particle. These experiments have also shown an 
asymmetry in the angular dependence of the coercivity. 



the coercivity reduction being more pronounced for ex- 
ternal fields along the particle axis than for fields directed 
perpendicularly to the sample. Both of these findings 
are in qualitative agreemento with the present theory. 
A quantitative comparison between theory and experi- 
ments is difficult for the presently available experimental 
data since the particles are irregularly shaped and often 
contain voids. Experiments on particles with a more per- 
fect morphology such as Cr02 or data of particles with 
various aspect ratios and diameters would clearly be de- 
sirable to further clarify the mechanism of thermally ac- 
tivated magnetization reversal. 

Let us now recall the various assumptions that have 
been made in the present theory: 

The cross-sectional area has been considered constant 
throughout the particle. This assumption leads to a con- 
tinuous degeneracy of the solution (f>s with respect to 
translations. In the case of a varying cross sectional area, 
the present treatment will still be approximately correct 
if the variations have a much shorter wavelength than 
the characteristic length scale S of the nucleus. How- 
ever, if the cross sectional area varies substantially, the 
saddle point energy will depend on the coordinate xq in 
( p.9[ ) and hence a whole class of energetically almost de- 
generate saddle points emerge. Such an extension of the 
present theory would predict that one single particle be- 
haves as if there would be a distribution of saddle point 
energies. Experimental resultsCJ indeed show deviations 
from an Arrhenius law involving a single energy barrier. 
They exhibit a decay of the magnetic moment of a single 
particle that is proportional to In t over several decades in 
t, a fact that is usually attributed to a distribution of en- 
ergy barriers. Such a behavior cannot be reconciled with 
the simple Neel-Brown picture which predicts a unique 
energy barrier for a single particle. 

In addition, we have focussed on nucleation in the inte- 
rior of the particle but we have neglected effects occurring 
at the particle ends: 

Since the nucleus (j)s describes a magnetization configu- 
ration merging asymptotically into the metastable state, 
the present nucleus may also be used to describe a sit- 
uation where the magnetization is pinned at the sample 
ends. In order for the present theory to hold, however, 
the pinning energy has to be sufficiently small that it can 
be overcome by the two domain walls propagating from 
the nucleation location to the sample ends. 

In the opposite case of free boundary conditions, i.e. 
M'(±L/2) = 0, there exists also the possibility that only 
one domain wall is nucleated at one sample end. This 
case can also be related to the present theory. In the ideal 
situation of a sample of constant cross section and an 
effective easy-axis anisotropy that extends to the sample 
end (at least within a distance smaller than the domain 
wall width), the saddle point structure 0^ restricted to 
the interval — 00 < a; < represents a domain wall which 
is nucleated at the sample end x — 0. Consequently the 
corresponding energy is half of the nucleus energy A£s- 

The theory as outlined in the appendix applies to the 
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regime of moderate to large friction. Since, however, the 
damping constant in magnetic systems is quite small, 
some estimates of the applicability range of the present 
theory are presented in the following. 



A. Validity of the Theory 

The principal existence of a lower limit of the damping 
constant for the present theory may be seen as follows. 
For a = 0, the linearized equations (|]^) do not describe 
the decay of the nucleus towards the stable state but 
rather a purely precessional motion which conserves the 
energy. Therefore, the corresponding decay frequency A+ 
is completely irrelevant for the nucleation rate for a — > 0. 

For very small values of the damping constant a, a 
completely different methodology would have to be ap- 
plied, since the nucleation does no longer correspond to 
a diffusion in configuration space but rather in energy 
space. Since the time evolution of the nucleus for ex- 
tremely spall a is expected to exhibit a "breathing" 
oscillationc3, a derivation of the corresponding Fokker- 
Planck equation would be an extremly difficult task. 
However, we shall see that the applicability range of the 
present theory extends to rather small values of a even 
for small nuclei and small cross sectional areas of the 
sample. Therefore we do not consider the underdamped 
theory any further. 

A criterion for the crossover between the present theory 
and energy diffusion has been given by Landauer and 
Swanson (Ref. |l^ and in Ref. ^ (see also Ref. ^)). The 
moderately damped theory may be applied if the energy 
loss during an (approximate) period of the motion near 
the saddle point exceeds ksT. 



Using the equations of motion (2.2) we obtain for the 
energy loss rate per area 



6£_ 



dx (9tm) 



.1) 



where m = M/Mq. Since (8.1 



) may also be expressed 
by spatial instead of temporal derivatives it is clear that 
the energy loss will be smallest for small nuclei. In 
this limit we may employ the spin wave approximation 
mj,, rriz <C 1, nix = — 1 -l- O(m^). The energy loss during 
one approximate period then takes the following form 



dt / dx {{dtmyf + {dtm,)^} (8.2) 



AA£ = -aA 



where ^ is the precession period. The r.h.s. of ( |8.2| ) 
is now evaluated approximately. First, we are only in- 
terested in leading order in a and thus we may use the 
conservative equations of motion for the evaluation of the 
integrand in (S.2). Second we discard any breathing ef- 
fects and neglect the exchange coupling of the magnetic 
moments such that the precession amplitude is given by 
the spatial distribution of the nucleus. Linearization of 



the equations of motion for a = 0, dtm = mx {d£o / dm) , 
with £o the energy (^]^) without exchange, then leads to 



dtniy = {1- h + Q ^)niz 
dtrriz = -(1 - h)my. 



(8.3) 



They describe an elliptical precession = m^sincjt, 
niy — —niy cos ujt with /niy — \/l — h/ y^l ~ h + Q^^ 
and uj = ^{1 — h){l — h + Q~^). Inserting this into 
( ^.2|) we obtain 

AA£ = -aA— [ rfa;(m°)2(l - h){l -h + ^). (8.4) 
LJ J 2 

Now, for small nuclei, we have h = sech^ i? ~ 1 — i?^. The 
precession is assumed to cover the nucleus structure and 



.0 _ 



2 sinh R cosh(x /S)/ (sinh^ R 



therefore we have 

cosh^{x/S)). Inserting this in ( ^.4[ ) and performing the 
integration the validity condition AA£ > kgT takes the 
following form in leading order R 



AS ^ d,TiQ-^l^'^(3AaR^ > 1. 



(8.5) 



For typical values of the constants as given above and 
Q^^ 0.2 we obtain the condition aR'^a 7 ■ 10^^ > 1 
where a is the cross section area measured in cm^. The 
limit is reached for e.g. a = 0.05, R = 0.4 {h = 0.84) at 
particle diameters 70A. Note that due to the neglection 
of the breathing contribution to the energy loss this rep- 
resents a lower limit for a and therefore the theory may 
be applied even for smaller values of a or samples smaller 
than indicated above. 
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APPENDIX 

In this appendix we present a derivation of the ex- 
pression for the n uclea tion rate starting from the Fokker- 
near the saddle point. 



Planck equation (3.8) near the saddle point. Since the 
dete rministic part of the dissipative linearized equations 
( |3.2[ ) constitutes a non Hermitian eigenvalue problem, 
there is no known set of eigenfunctions to expand in and 
therefore traditional formulations which are based on an 
expansion into mode amplitudes cannot be applied im- 
mediately. In the sequel we shall show, however, that 
the within a functional formulation, the derivation of the 
nucleation rate can be. carried out in close analogy to the 
methods of Kramerslij and its extension to many degrees 
of freedom by LangerQ'^'') . In the limit of large friction the 
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result reduces to the theorieS of Brinkman, Landauer- 
Swanson and Langer. 

The IpfliSic principles of the method go back to 
Krameralj. We calculate the stationary flux of a nonequi- 
librium distribution across a surface transverse to the 
unstable direction at the saddle point. A main differ- 
ence to finite dimensional problems is the existence of 
the Goldstone mode with zero energy which reflects the 
continuous degeneracy of the nucleus (fis with respect to 
rigid translations. 

Our goal is the construction of a stationary nonequi- 
librium probability density obeying the boundary condi- 
tions g 1 near the metastable state and £i ~ beyond 
the saddle point. To this end we factorize the desired 
probability distribution as follows: 



g = QcqF. 



(Al) 



The key assumption is now to let be a function of one 
coordinate u only. F is such that g is a normalizablc func- 
tion. In the vicinity of the saddle point, the coordinate 
u is a linear functional in ^ 



= J dxJ2Ujix)^j{x), 



(A2) 



with {Ui,U2)^ {Uf,UP). After insertion of (|Al|) with 
([A^) into (3.S) and using the stationarity of g we obtain 



-'Hi'4>i{x)Uj{x) 



dF 



du 
= 0. 



(A3) 



In order for (A3) to become a proper differential equation 
in u alone wc must have 



dxY^ M,jH^^,{x)Uj{x) 



KU, 



where 



7, K 



dx^A'hjUi{x)Uj{x) = 



7K, 



(A4) 



(A5) 



are real constants. Since^ h'om (3.5), 

— ctJ2iUf, the condition (A£) amounts 



t o a normalization o f th e Ui and we have 7K > 0. Using 
( A4) and (A5), Eq (|A^ reduces to the same differential 
equation as in the casetj of one degree of freedom: 



dF 
du 



■7 



d^F 
dv? 



0. 



(A6) 



which is integrated with the boundary conditions 
F(— 00) — 1, -F'(oo) = 0, i.e. the vicinity of the 
metastable state [u — > —00) is characterized by thermal 
equilibrium while the probability distribution vanishes 



beyond the sadlle point {u 00). Since g has to be nor- 
malizablc we must have 7 < and hence k < 0. -F is 
then given by 



F{u) 



1 r , / 

— , / du exp < — -- — r 
y/2^\Ju I 2|7| 



(A7) 



Inserting (A7) and ( |A5| ) into (3^) we obtain for the cur- 
rent near the saddle point: 



1^ 
2|7l 



exp <! \ goq M^jUji 



X 



(A8) 



We now return to the evaluation of k. Equation (A4) is 



fulfilled if 



(A9) 



Note that this implies that U'^{x) has no component 
along the zero frequency mode since [Xi' ^U^) 



0. This is physically plausible 



since the dynamical instability of (f>s is only associated 
with a shrinking or expansion of the nucleus but not with 
a pure translation of the nucleus as described by Xi'' (x) ■ 
Because is a positive operator there is no restriction 
on the functions U'p{x). Therefore (A9) can be inverted 
to give 



(AlO) 



where (TY*'^)"-'^ acts only on the subspace {/|(/, Xi"^) = 
0}. Eqn ( Alpp and hence ( A9) are solved by putting Ui = 
Tiitlif where ij)^ = {ip+,p+) is the dynamical unstable 



mode obeying (3.2). Therefore we have 



= -A+ < 0. 



(All) 



The escape rate F is now obtained by integrating the 
flux transverse to the unstable direction, e.g. over the 
manifold u = 0: 

T = [ VifVp [ dxJ2Ut{x)Ji{x) 




(A12) 



where we have used ([A8|) and A+ is determined by the 
eigenvalue equation (3.2). The constant 7 will cancel 
in the final result as we shall see below, goq has to 
be evaluated on the hyperplane u = in the vicinity 
of the saddle point. According to ( 3.10 ) it is given by 



Pcq ~ exp{— /3y^£s }/Z with . Z is determined by the 
normalization of gcq in the vicinity of the metastable 
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state. Upon insertion of (3.1C) and usin g the Fourier 
representation of the (5-function, eq. ( A12 ) reads 



r = A. 



exp |-/3^[^ J dxipH^^if dxpW^p] 



(A13) 



For the evaluation of the functional integrals we use the 
self adjointness of TL^'^ and TC^^ acting on functions with 
periodic boundary conditions on [— -j, ■§■]• Therefore they 
have an orthogonal and complete set of eigenfunctions 

x7.x7 with n'^x7 = e7x7, n^'xtF = E^^TxtF, 

where i', ^ denote the bound states v, fi — 0,1,(2) and 
scattering states v, fi = k as well. Therefore we can ex- 
pand 



(p{x) 
p{x) 



Vi.x7{^). 



(A14) 
(A15) 



where X-fe = Xk-> ^'^^ '4^viVv are complex expansion co- 
efficients. Since ip and p are real we have Lp^ = i^al^ 
and pv = P*-i,- Expanding similarly u = 'Yli^'i * '^v ^ 
Yl, UJl *Pt^ we obtain 



r 




Wdip,, jWdp 



X exp 



exp <^ iqll^ipo 



PA 



E'o'^lvl 



(A18) 



The expression det H stands for the product of all eigen- 
values of Ti.. The double prime on the determinant and 
the summation indicates that the terms v = 0,1 corre- 
sponding to the unstable mode and the zero mode are 
omitted. 

We now turn to the evaluation of the zero mode. 
To this end we remember that a pure translation of 
the nucleus can be described by (/)s{x + dx) — 4>s{x) — 
(j)'g{x)dx or since x'l' oc dcp^/dx we can equally well write 
(j3s(x-\-dx) — 4>s(x) = xVi^) dfi- This allows us to replace 
the integration over the zero mode by an integration over 
X. Since Xi'^ is normalized to unity and reinstating the 
integration measure (A17) of {dipi) we have 



{dipi) = V — 




where 



(A19) 



(A20) 



Here we have made use of (2.15). After performing the 
q- integration ( [AlSj ) we have to convince ourselves that 
the remaining integral converges. This is e nsure d by the 
following relation, which follows from (A£) and ( AlC ) 



^ y" dxU^{x)'H'^^U,{x) = -yi\(iA < 0. (A21) 



X exp 



PA 



EfFp:p,) 



(A16) 



where the prime on the sum indicates that the integrand 
is independent of ipi, since E7 = Uf = 0. To simplify 
notation we n ow choose integration measures (dip^) and 
(dp^) in (A16) and in Z such that e.g. 



(dpo) exp 



PA 



^0 Pf) 



{dpl){dpk)eM-0AElPplpk} 




(A17) 



But note that we have to restore -y/ PA/2tt in the in- 
tegration over the zero mode which does not contain a 
Gaussian and note also that the integration measure in 
q is the usual one. We now perform the integrations in 
if I, and Pf^ except for the amplitudes of the zero mode ipi 
and that of the unstable mode, ipo, to obtain: 



A, 



7l e 



1 



1 



2tt Z 



dq 



X / (d^o) I ^ exp 



Vdet" VdetWp 

2 „ 



{dipi 



2PA^^'^ E, 




Using the expansions of Ui in terms of the x7^'^ and 
recalling that [/f = we see that the exponent is in 
fact negative. We are now left with the evaluation of Z. 
Since g^q is strongly peaked at the metastable state we 
can perform a Gaussian approximation 



Z^ J VLpVpexp^^-PA[^ J dxipT-C- 



/ dxpTT'Pp] 



(A22) 



where in contrast to the above (p and p now describe 
fluctuations out of the metastable state. They are defined 
hy (j) = TT + ip and 6 = 7r/2 — p, where \(p\, \p\ <C 1. This 
integral is performed analogously to the previous one: (p 
and p may be expanded into the (plane wave) eigenstates 
oiTC^'f ,T-C^P . The i ntegr ations then are all Gaussian and 
using the measure (A17) for the integrations, we finally 
obtain 



1 



1 



Vdet H"'^ Vdet Ti^P 



(A23) 



After having performed the g-integration we can now 
carry out the fina l Ga us sian integration over the unstable 
mode ipQ. Using ([A5|), (A18) leads to the final result: 
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-PAS, 



det WP 



(A24) 



where det' | | denotes the product of the modulus of 
the eigenvalues with omission of the zero mode. Apart 
from the dynamical prefactor the result looks as if we had 
evaluated in Gaussian approximation the ratio Z^/Zs of 
the partition functions at the metastable state and the 
saddle point, respectively, with the unstable mode ren- 
dered to a stable one. Or loosely speaking, the nucleation 
rate is proportional to the imaginary part of the ratio 
Zm/Zg- For later reference, we express the result ( A24 ) 
in terms of the eigenvalues El'f ,E^P: 



r 




1 



27r ,/\E^^\E^^EfEf 



(A25) 



where according to ( A20| ) C = L^/El {£s is given by ( 2.15| ) 
) and A+ is determined by (3^). Note that in our case 
of two equivalent saddle points 0^ we have to multiply 
this final result by a factor of two. 
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FIG. 1. Various anisotropy configurations axes that can 



be described by the energy density (2.5). The sample cross 
sectional areas are assumed to be sufficiently small such that 
transversal variations in the magnetization are suppressed. 



FIG. 2. The spatial variation of the nucleus is shown for 
i) fields close to the anisotropy field and ii) small fields. In 
these pictures, the magnetic chain is taken along the easy-axis. 



However, the model (2.f:) equivalently applies to all situations 
shown in Fig ^ 



FIG. 3. The reduced prefactor in the overdamped limit is 
shown as a function of the parameter R for different values 
of the hard-axis aniso trop y. Th e lin es ( — ) and ( ) are the 



asymptotic formulas (5.5) and (5.3) respectively. 



FIG. 4. The decay frequency of the nucleus is shown for 
different values of R and a. The dots are results of a numer- 



formula (6.2) 



ical solut ion of (3.2) and the solid line is the approximation 



FIG. 5. Numerical results for the reduced prefactor (4.55) 
for moderate damping as a function of R. The dashed lines 
represent the asymptotic formula (^.Tj). 



FIG. 6. The total nucleation rate F is shown as a function 
of the reduced external field for various particle diameters. 
The material parameters are chosen as in Sec. VIII. The 
dashed curves indicate the results of the Neel-Brown theory. 
The particle aspect ratio is assumed to be 15:1. 
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TABLE I. Summary of the results for the prefactor fl (4.56) in the moderately damped limit for i) large effective hard-axis 
anisotrop y Q ~^5 3> 1, (5 = cothi?), and ii) small fields i? ^ 1 (/i = sech^i?). Note that the underlying dimensionless units are 
given by (|2~ 



A4 



sifi b 



'S' > l",i?< 1 



r'<5' > 1, 

/QR/a > 1'' 



-'s'^ > r,i?< 1, 



i^tanh^/^i^sinhii 



+ 



^(f)'[Q-2-2Q-i£on 



y/3Q-^R 



-1 fP^V 



(Enl) 



3(a + 



(8.9) 



(6.9) 



(6.5) 



a) i? > 1 



i?> 1, 



3^ tanh^/^ i?sinhi? 



(6.7) 



(q + a"')e " 



-1 



8(a + a-^)e- 



(5.2) 



(6.7) 



(6.2)^ 



(4.2e) 



In this limit, det-W'"f/dcfW' = 1, cf. ( [4.3lD 
''has to be evalu ated numerically 
Compared to (6.2), small terms of the order C'((Q_Eq''')^) have been dropped. 



typically realized in experiments 

'For a 00, these results merge into (5.f:), (5.3) obtained in the overdamped limit 
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